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Abstract 

We study, both analytically and by numerical modeling the equilibrium prob- 
ability density function for an non - linear Levy oscillator with the Levy index 
a, 1 < a < 2, and the potential energy x 4 . In particular, we show that the 
equilibrium PDF is bimodal and has power law asymptotics with the exponent 
-(a + 3). 
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1 Starting equations 

Recently, kinetic equations with fractional derivatives have attracted attention as a 
possible tool for the description of anomalous diffusion and relaxation phenomena, see, 
e.g., recent review [ij, semi - review papers [0,0,13 and references on earlier stud- 
ies therein. It was also recognized |Q, || that the fractional kinetic equations may 
be viewed as "hydrodynamic" (that is, long - time and long - space) limits of the 
CTRW (Continuous Time Random Walk) scheme J7|, which was successfully applied 
to the description of anomalous diffusion phenomena in many areas, e.g., turbulence 



||, disordered medium ||, intermittent chaotic systems |TD[, etc. However, the kinetic 
equations have two advantages over a random walk approach: firstly, they allows one to 
explore various boundary conditions (e.g., reflecting and/or absorbing) and, secondly, 
to study diffusion and/or relaxation phenomena in external fields, both possibilities 
are difficult to realize in the framework of CTRW (we point, however, to the paper 



TTfl , in which a fractional kinetic equation was obtained from generalized CTRW). 
Fractional kinetic equations can be divided into three classes: the first one, describing 
Markovian processes, contains equations with fractional space or velocity derivatives 
and the first time derivative, the second one, describing non - Markovian processes, 
contains equations with fractional time derivative, and the third class, naturally, con- 
tains both fractional space and time derivatives, as well. In this paper we deal with a 
one - dimensional kinetic equation belonging to the first class, namely, with the Frac- 
tional Symmetric Einstein - Smoluchowski Equation (FSESE), which, from one hand, 
is a natural generalization of the diffusion - like equation with the symmetric fractional 



space derivative [H],|jl2| and, from the other hand, is a Markovian generalization of 



the Einstein - Smoluchowski kinetic equation, which describes a motion of a particle 
subjected to a white Gaussian noise in a strong friction limit, see, e.g., From this 



point of view, the FSESE describes a motion of a particle subjected to a white Levy 



noise, also in a strong friction limit JI4 



In dimensionless units the one - dimensional FSESE has the form 
d f d d a f 

~K7 — ( F f) + TTT^ , t > 0, X eR, f( X , 0) = 5{X) , (1.1) 

at ox o \x\ 

where f{x,t) is the probability density function, F is the deterministic external force, 
a is the Levy index, < a < 2, d a /d \x\ a is the symmetric fractional space derivative 
[2, 11], which is defined, for a "sufficiently well - behaved function <fi{x) iL through its 
Fourier - transform <p{k) as 

'(*)-*-- |*|° (1.2) 



d \x\ 

or in terms of the Riemann - Liouville derivatives as 



d\x\ a 2cos(7ra/2) 
where a > 0, a ^ 1, 3, ... , 

D<*d>(x) = - ±lWW_ D «a ) (x) = - - 1 ^ X ' )dX ' (14) 

+nX) r(l - a) dx J (x- x'Y ' ^ ] T(l-a)dxJ (x' - z)° 1 ' 

— OO X 

for < a < 1. For a > 1 



±rv ; Tn-a)dx n 1 



o 



n — 1 < a < n E N. The derivatives (1.4) and (1.5) are characterized in their Fourier 
representation by 

D±(f)(x) (^ik) a (f)(k), (1.6) 



where {^ik) a = \k\ a exp ^^-sgnk^j . One can easily recover Eq. (1.2) by combining 
Eq. (1.3) with Eq. (1.6). A detailed theory of the Riemann - Liouville and other forms 
of fractional derivatives is presented in ||15|| . 

The solution of Eq. (1.1) in the force - free case, F = 0, is known to be a probability 
density function (PDF) for an a— stable symmetric process. The solution of Eq. (1.1) 
was also obtained for the constant field, F = const on the whole axis, — oo < x < oo, 
and for a linear Levy oscillator [16|]. As far as the authors know, Eq. (1.1) was not 



considered for a more complicated potential fields. In the present paper we study the 
equilibrium solution of the FSESE for non - linear Levy oscillator with the potential 
energy 

r 4 dU 

The models with the potential energy (1.7) play an important role in the theory of 



dynamical chaos [17| and in the theory of Brownian motion in an open auto - oscillation 



systems |18|| , and have various applications One may expect that the models of non 
- linear Levy oscillator will also possess an important place in the theory of the systems 
influenced by non - Gaussian Levy noises and obeying fractional kinetic equations. For 
our purposes we pass to the equation for the characteristic function f{k) , 

oo 

f(k) — dxf{x)exp(ikx). (1.8) 



For the equilibrium state f(k) obeys the equation, which follows from Eqs. (1.1) and 
(1.7): 

^ = sgn(k)\kr 1 f eq (k), (1.9) 

where the index "eq." denotes equilibrium solution. Below in the paper we omit 
subscript "eq" for brevity. The characteristic function obeys the following conditions. 

1. f (k) = f* (k) = f (—k) ,where the asterisk implies complex conjugate. The first 
equality is a consequence of the Khintchine theorem about reality of the characteristic 
function for the symmetric PDF, whereas the second equality is the consequence of the 
Bochner - Khintchine theorem about positive definiteness of the characteristic function. 
In fact, one can see that the solution of Eq. (1.9) obeys the condition 1. 

2. The boundary conditions for f(k) are as follows: 

/(0) = 1, /(±oo) = 0, ^ = 0- (1.10) 

oo 

The first property stems from normalization of the PDF, J f(x)dx = 1. 

— oo 

The second property stems from the existence of the PDF. 

As to the third property, we remind that the integer moments of the PDF (if exist) 
are connected with the derivatives of the characteristic function at k = as 

, „. 1 (0) 

(x p ) = - „ , P= 1,2... 

N ' %p dkP 

The third property is the consequence of this theorem: 

df ip) (0) n . , , 

——— = 0, p=l, 3, 5, 



because the PDF is a symmetric function, and hence, all odd moments are equal zero. 
The last equality is valid for those odd p, for which the p - th moments of the PDF 
exist. 

To make results more transparent, the rest part of the paper will be organized 
as follows. In Sect. 2 we recall the results for non - linear Brownian oscillator, a = 
2. In Sect. 3 we get an equilibrium solution for an non - linear Cauchy oscillator, 
a — 1. In Sect. 4, which is the main one, we present analytical and numerical results 
on equilibrium solutions for non - linear Levy oscillators with the Levy indexes lying 
between 1 and 2, Brownian and Cauchy oscillators being their particular cases. Finally, 
in Sect. 5 we present the discussion and summarize the results. Some subsidiary studies, 
which support the results obtained in Sec. 4, are presented in Appendix A. An integral 
which is necessary for estimating asymptotics of the PDF is calculated in Appendix B. 



2 Equilibrium solution for a non - linear Brownian 
oscillator 

In this Section, which is presented here mainly for the methodical purposes, we remind 
the equilibrium solution for non - linear Brownian oscillator, a = 2. It is well - known, 
that in this case there is no need to pass to the equation (1.9); on the contrary, the 
starting point is the stationary equation, which follows from Eqs. (1.1) and (1.7): 

*^ --§'<•>■ <"> 

It has the Boltzmann solution 



f(x) = Cexp(-U(x)), (2.2) 

where the constant C is determined from normalization condition. For a non - linear 
quartic Brownian oscillator 

c =rm- (2 - 3) 

We are also in position to get the characteristic function by making Fourier trans- 
form of equilibrium solution (2.2), (2.3) and, then expanding expiikx) into the power 
series and integrating over x each term separately: 

oo 

/>)=j>A* (2.4) 

j=0 

where 

(-i)^rm 
°>- m r(i/4)- (2 ' 5j 

In Sect.. 4 we show that this solution is a particular case of that for a non - linear 
Levy oscillator. 



3 Equilibrium solution for non - linear Cauchy os- 
cillator 



For this case we start from Eq. (1.9) with a — 1. The solution is 
.This solution can be expanded into the power series of \k\ as 

oo 

/(*)=!> |*|'\ (3-2) 

j=0 



where 

1 



2j! 



j \ (2-ni \ ( i \ /4m 
-7= ex P -^-J + 1 + ex P -Z-j 



(3.3) 



V3; v 3 j j v v/s; v 3 

Making an inverse Fourier transform of Eq. (3.1), we get equilibrium PDF for the 
Cauchy oscillator, 

f( X ) = — \ jr. (3.4) 

J W 7T(1 - X 2 + X 4 ) V 7 

oo 

Clearly, /(x) > and f f{x)dx — 1, so / is a PDF. 

— oo 

Equation (3.4) clearly indicates strongly non - Boltzmann character of the equilib- 
rium PDF for non - linear Cauchy oscillator. In Fig.l both PDFs for Brownian and 
Cauchy oscillators are shown in a linear (at the top) and in a semi - logarithmic (at the 
bottom) scales. The two important distinctive features of the equilibrium PDF for the 
Cauchy oscillator are (i) the bimodality, and (ii) power law asymptotic at x — > ±oo. 
The latter property is clearly visualized in a linear scale, whereas the former is better 
shown in a semi - logarithmic scale. It appears that both features are inherent not only 
to a non - linear Cauchy oscillator, but also to equilibrium PDFs of non - linear Levy 
oscillators with the Levy indexes a , such that 1 < a < 2. In the next Section we get 
the solution for an arbitrary a between 1 and 2 and clearly demonstrate this fact. 



4 Equilibrium solution for non - linear Levy oscil- 
lator, 1 < a < 2 

We turn to the solution of Eq. (1.9) with the boundary conditions (1.10). One can 
convince himself that the solution of Eq. (1.9) with the first and the third conditions 
from Eqs. (1.10) being taken into account can be represented as 

/>) = S 1 + aA; 2 S 2 , (4.1) 

where 

oo 

E^l + E^l^, (4-2) 



£ 2 = i + £W (a+2) . (4.3) 

3=1 

and the coefficients aj, bj are obtained by inserting Eqs. (4.2) and (4.3) into Eq. (1.9) 
and equating the terms of the same powers of k in the right - and left - hand sides: 

a j3 (a + 2)(ja + 2j - l)(ja + 2j - 2) = (4.4) 

bjj(a + 2){ja + 2j + l)(ja + 2j + 2) = b^, (4.5) 
j > 1, a = b = 1. 

Obviously, the CLj cLS well as bj tend to zero extremely fast. 

The terms Ei and k 2 T,2 are, in fact, two independent particular solutions of Eq. 
(1.9). Since the condition at the infinity from Eq. (1.10) has not been employed yet, 
the general solution thus depends on an arbitrary constant a. We define it numerically 
by demanding 

f(k - oo) - 0, (4.6) 

that is, 

Our numerical simulations show that with k increasing the value of a rapidly reaches 
the constant value, which, of course, depends on a. Further increase of k allows us to get 
a with higher accuracy, that is, with more significant digits. It is also worthwhile to note 
that the radius R of convergence for both power sets in numerator and denominator 
of Eq. (4.7) is infinite. This fact can be easily shown with the help of the Cauchy - 



Hadamard theorem (see, e.g., IpOfl , p. 300), according to which, e.g., for Ei 

— = lim \aj\ lf \ 

It J^co 

where the bar denotes the largest limit for the sequence j|aj| 1//jf 

In Fig. 2 the obtained solution for f(k) is shown in a linear (at the top) and in a 
semi - logarithmic (in the bottom) scales, a = 1.7. One can see an oscillatory character 
of the solution at large k; this property can not be visualized directly from the power 
series expansion. Therefore, we also get a large k— asymptotic of the solution to Eq. 
(1.9). The derivation of it is presented in Appendix A in detail. There we demonstrate 
that (i) the asymptotics has an oscillatory character and an exponentially decreasing 
amplitude, and (ii) the period of oscillations coincide with high accuracy with that 
obtained numerically from power series expansion, see Fig.2. Thus, with the help of 
Appendix A we have an independent evidence of the correctness of our procedure and 
of the solution presented in this Section above. 

Another evidence of the correctness of the obtained results stems from the com- 
parison of the general expressions (4.1) - (4.5) with those for the two particular cases 
a = 2 and 1, presented in Sect.. 2, see Eqs. (2.4), (2.5), and in Sect.. 3, see Eqs. 
(3.2), (3.3), respectively. Indeed, let us consider Brownian oscillator at first. By com- 
paring expansion (4.1) at = 2 with that given by Eq. (2.4) one can see that the odd 



coefficients C2j and C2j+2 given by Eq. (2.5) have the same recurrent relation as the 
coefficients dj and a J+ i given by Eq. (4.4) have, whereas the even coefficients C2j+i 
and C2j+3 given by Eq. (2.5) have the same recurrent relation as the coefficients bj and 
bj + i given by Eq. (4.5) have. Now, let us compare general expressions with those for 
the Cauchy oscillator. By comparing expansion (4.1) at a = 1 with that given by Eq. 
(3.2) one can see that the coefficients c^j and c 3 j +3 in Eq. (3.3) obey recurrent relation 
between the coefficients aj and Eq. (4.4), whereas the coefficients c 3j _i and 

%+2 obey recurrent relation for the coefficients bj and b j+ i, see Eq. (4.5). At last, 
C3 J+ i = 0, j > 1, which is also in agreement with general expansion, if we set a = 1 
in Eq. (4.1). Thus, we may conclude that the general expansion (4.1) - (4.5) is in 
agreement with the particular cases of Brownian and Cauchy oscillators. However, we 
have to take in mind that for the Brownian and Cauchy cases we are able to evaluate 
coefficient a analytically. On the contrary, for an arbitrary we have to estimate this 
coefficient numerically with the use of Eq. (4.7). Of course, the results of numerical 
and analytical estimates coincide for a = 1 and 2. 

Now we consider two important properties which have been already discussed for 
the particular case of the Cauchy oscillator in Sect.. 3, namely, power law tails and 
bimodality. Consider power law tails at x — > ±oo at first. These asymptotics are 
determined by the first non - analytical term in the power series expansion (4.1), that 
is, the term ai|/c| Q+2 . By making an inverse Fourier transform of this term, we get 

oo 

f(x)ttai / — \k\ a+2 exp(— ikx) , x — > ±oo, (4.8) 
J 2-k 

— oo 

where a\ = [a(at + l)(a + 2)] _1 . This integral is calculated by passing to the 
complex plane, the detailed derivation is presented in Appendix B. The result is 

„. . sin(7ro;/2)r(a) , , , 

/(*) * | |o+3 ,*-±°o. (4.9) 

7T \X\ 

It follows from Eq. (4.8) that the equilibrium PDF has a power law tail, f(x) oc 
j^j-O+s) ^ an( ^ fays the integer moments of the order greater than 3 diverge. This 
behavior is strikingly different from that of a non - linear Brownian oscillator. The 
"long tails" can be explained qualitatively, if we turn to the Langevin description of 
the Levy oscillator, The Langevin approach relevant to the FSESE [13] implies that 
the non - linear overdamped oscillator is influenced by "white Levy noise" £(£), whose 
PDF behaves as |£| -1 ~"at |£| — > oo . These "long tails" imply that the large absolute 
values of the noise occur frequently, which, in turn, lead to large increments of the 
coordinate. However, it is also clear that the PDF of the coordinate x must fall off 
more rapidly at x — > oo than the PDF of the noise £, because of the presence of the 
potential well, which prevents x from "escaping" rather far from the origin. 

In Fig. 3 equilibrium PDF is shown by solid lines in a linear (at the top) and 
semi - logarithmic (at the bottom) scales. The PDF is obtained by an inverse Fourier 
transform of the characteristic function shown in Fig. 2, a = 1.5. The dashed lines 
indicate asymptotic (4.9). One can see, especially from the semi - logarithmic plot, that 
the asymptotics is a good approximation beginning from k equal nearly 2. In this figure 
the second important property, namely, bimodality is clearly seen in a linear scale. In 
Fig. 4 the profiles of equilibrium PDFs (obtained by an inverse Fourier transform of 
corresponding power series in k -space) are shown for the different Levy indexes from 



a = 1 at the top of the figure till a = 2 at the bottom. It is seen that the bimodality is 
most strongly expressed for a = 1. With the Levy index increasing the bimodal profile 
"smooths" and, finally, it turns to a unimodal one at a = 2, that is, for the Boltzmann 
distribution. 

We also perform numerical simulation based on numerical solution of the Langevin 
equation 

f = - 3+ «- < 4 - 10 > 

where £(t) is a white Gaussian noise or a white Levy noise, whose generator is described 
in detail in our previous papers devoted to the studies of self - affine properties of 
ordinary and fractional Levy motions [2~i~l,|[22||. A discussion on equivalency of the 



description of a stochastic system with the help of the Langevin equation (21) and 



fractional kinetic equation (1) is presented in Ref. | Hfl ; for the Brownian motion this 
problem is discussed, e.g., in Ref. fll3"l , in detail. 



In Fig. 5 the results of numerical modelling are presented for the Brownian oscillator, 
a = 2 (above) and the Levy oscillator, a — 1.1 (below). At the left the trajectories x(N) 
are presented , where N = t/At is the number of time steps in numerical modeling, t is 
the length of a single step, At << 1. The result of numerical solution of the Langevin 
equation must not depend on At; for the Brownian oscillator this requirement is fulfilled 
at At < 10~ 2 , whereas for the Levy oscillator At < 10 -3 . We also studied the time - 
dependence of the second moments and fix when the moments become constant, thus 
indicating equilibrium state. From the left figures a clear difference between trajectories 
of Brownian (above) and Levy (below) oscillators is seen: there are large "jumps" on the 
figure below, which are due to existence of large "pushes" from an external Levy noise 
£ or, equivalently, due to power law asymptotic of the PDF of the Levy noise £. Now we 
turn to the right figures. The designations are as follows. Thick solid line 1 and dotted 
line 2 show the potential well and its curvature, respectively, in conventional units. Thin 
solid line 3 shows the Boltzmann distribution (2.2). The black points depict the PDFs 
obtained in numerical simulations by statistical averaging over 50 trajectories each of 
20.000 steps. Finally, the power law asymptotic (4.9) are depicted by a thin solid line 4. 
It is seen from the figure above that the PDF obtained in numerical simulations agrees 
quantitatively with the Boltzmann PDF, whereas the figure below demonstrates drastic 
difference between the Boltzmann PDF and numerical PDF. The latter has long power 
asymptotic, which start far away from the maximum of the potential well curvature, 
and whose exponent is close to that obtained theoretically. This conclusion is confirmed 
by the results of numerical modeling for the Levy oscillators with different Levy indexes, 
see Fig. 6. In this figure the black points depict the values 7 of the exponents of the 
asymptotic |a;|~ 7 - of the equilibrium PDFs for a = 1.1, 1.3, 1.5, 1.7 and 1.9. The values 
7 of are estimated as a tangent of a slope angle of a rectilinear asymptotic in a double 
logarithmic scale. The statistical averaging is over 50 trajectories, each consisting of 
20.000 time steps, At = 10 -3 . Dotted line depicts theoretical dependence 7 = a + 3 of 
the asymptotic exponent versus a. The results of numerical simulations coincide with 
analytical estimates within the error limits. 

In Fig. 7 the results of numerical modeling are presented in more detail for small 
values of x, when power law asymptotic is inadequate. Similarly to Fig. 5, thick solid 
line 1 and dotted line 2 show the potential well and its curvature, respectively. Solid 
line 3 indicates PDF obtained by an inverse Fourier transform of the power series in 
k— space. The black points depict PDF obtained in numerical simulations by statistical 



averaging over 20 trajectories , each of 10~ 5 steps. The Levy index is 1.2. It is clearly 
seen that both PDFs are in qualitative agreement and the both are bimodal. Therefore, 
from Figs. 5 - 7 we may conclude that the results of simulations based on numerical 
solution of the Langevin equation confirm qualitatively and even quantitatively the 
results based on analytical solution of Eq. (1.9) for the characteristic function. 

5 Discussion and Results 

In this paper we study, both analytically and numerically, the properties of equilibrium 
PDF of a non - linear (x 4 ) Levy oscillator, that is, the oscillator which is subjected 
to a white Levy noise obeying a Levy stable probability law. We restrict ourselves 
to the case of the Levy indexes such that 1 < a < 2. It is known that the Levy 
stable distributions (as the Gaussian one, which corresponds to a equal 2) appear in 
problems, whose result is determined by the sum of a great number of independent 
identical factors. Since the Brownian oscillator is subjected to a white Gaussian noise, 
the Levy oscillator is a natural generalization of a Brownian one. 

For the analytical studies the starting equation is the so - called fractional symmetric 
Einstein - Smoluchowski equation, which contains fractional symmetric space derivative 
and is a natural generalization of the kinetic Einstein - Smoluchowski equation used in 
the theory of Brownian motion. 

The main results are as follows: 

1. We get analytically the characteristic function of the equilibrium PDF in the 
form of a power series. Its inverse Fourier transform, realized numerically, allows us to 
obtain the PDF. The two main distinctive features of the equilibrium PDFs for non 
- linear Levy oscillator with a ^ 2 are (i) the power law asymptotic at large x, and 
(ii) bimodality. Both features imply that the PDFs for non - linear Levy oscillators 
are strikingly different from Boltzmann distribution, which is equilibrium PDF for the 
Brownian oscillator. 

2. The power law asymptotics is determined by the first non - analytical term in 
a power series expansion of the characteristic function. We find that the asymptotics 
behaves as \x\~( a+3 h 

4. The bimodal profile of the equilibrium PDF is clearly visualized after an inverse 
Fourier transform of the power series expansion for the characteristic function. Two 
maxima are most strongly expressed for a = 1 (Cauchy oscillator). With the Levy 
index increasing the profile "smooths" and at a = 2 (Brownian oscillator) turns to 
profile with a single maximum. 

5. We make three independent verifications of the solution for the characteristic 
function. 

First: we find that the power series expansion for the characteristic function of the 
Brownian oscillator (a = 2) is in agreement with the solution for 1 < a < 2, if we set 
a = 2 in the latter. 

Second: the particular case of the Cauchy oscillator ( a = 1) admits complete 
analytical study, which allows us also to check the agreement between power series and 
to demonstrate power law asymptotics and bimodality of the PDF. 

Third: we obtain the solution of the equation for the characteristic function at large 
values of the argument k and show that the periods of oscillations of this asymptotic 
solution in A;— space coincide with high accuracy with those obtained from the power 
series. 



Therefore, all three particular verifications testify to correctness of our approach in 
the general case, 1 < a < 2. 

6. We also perform numerical simulation based on numerical solution of the Langevin 
equation for a non - linear oscillator subjected to a white Gaussian noise. The equilib- 
rium PDFs obtained in simulation show quantitative agreement with PDFs obtained 
analytically. 



At the end we note that the general case of a non - linear Levy oscillator (x 



2n 



71 



1, 2, ...) can be treated in a similar way. In particular, it can be shown that in this case 
the power law asymptotics have an exponent — (a + 2n — 1). 
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APPENDIX A. Asymptotics of equilibrium characteristic func- 
tion at large k 



We start from Eq. (1.9) for k > 0. Using the transformation |23 



V(0 = k b/3 f(k), i = ck 1+b '\ (A.l) 
where b = a — 1, and c is an arbitrary positive (for definiteness) parameter, we get 



^ " 

t 3 v 



(1 - u 2 )^ 
At large Eq. (A. 2) reduces to 



av 



1 - ~r? v = 



a + 2 



(A.2) 



whose solution, which tends to zero at large is 



V(0 = Ci exp 



c 



1 iVS 
1 + ^ 



+ C 2 exp 



c 



1 
2 



(A3) 



where C\ and Ci are arbitrary constants. Returning to f(k) , we get at large k 



f(k) = cat (q - 1)/3 



exp 



2(a + 2) 



fe(a+2)/3 



X COS 



3V3 
2(a + 2) 



fc (a+2)/3 _ q 



(A.4) 



It depends on two unknown constants, C and 8, since we use the condition at 
infinity only. We can check this formula for two particular cases, a = 1 and 2. 

1. Cauchy oscillator, a = 1. 

By comparing Eqs. (A.4) and (3.1) we get C = 2/y/S , 9 = %/6 . 

2. Brownian oscillator, a = 2. 

We explore asymptotics of the following integral at large k, which can be obtained 
with the help of the saddle - point method, see, e.g., pi| : 



CO 

J dxm-p(-x 2m +ikx) 



^Jm(2m — 1) 



( 2m )(™-l)/(2m-l) £.-(m-l)/(2m-l) x 



exp ( k 2m /^c m cos^^-] cos ( fc 2 -/(2-D in _ ^rn-i; 

P V 2m- 1 J V 2m -1 2 2m - 1 

(A.5) 



where m — 1,2, 



2m V ; 



With the help of Eq. (A.5) we get, also using Eqs. (2.2), (2.3), 

oo oo 



X 4 



4 



2 V^ l7 .|-l/3_ f 3 4/3 \ / 3^ 4/3 



fcp^exp ~ |A;| 4/d cos ^ |A;| 4/d - - (A.6) 



>/3r(l/4) \ & J ^ 8 1 6 

at | A; | >> 1. By comparing Eqs. (A.6) and (A.4) we get C = 2^/(^(1/4)) , = 
tt/6. 

We also make a comparison between the period of the asymptotics (A.4) and the 
period of oscillations of the solution (4.1) - (4.7). Firstly, we are able to get with high 
accuracy the values kj, j — 1, 2, 3, 4, 5, at which the solution given by Eqs. (4.1) - (4.7) 
is equal zero, see also Fig. 2. Then, we insert kj into the cosine in Eq. (A.4) and 
estimate 

o fo l.(«+2)/3 l>+2)/3 

1 " J^,— A = S j+lj . (A.7) 

2(a + 2) 7T 3+ ' 3 y ' 

It is seen, that Sj+ij can serve as a measure of difference between zeros of f(k) 
estimated from Eqs. (4.1) - (4.7) and those estimated from Eq. (A.4). For example, 
for the characteristic function with a = 1.7, which is shown in Fig. 2, we get #2,1 = 
56 • 10~ 4 , S 3:2 = 17 • 10" 4 , <5 4 , 3 = 8 • 10~ 4 , <5 5)4 = 5 • 10" 4 . These results demonstrate 
the smallness of Sj+ij and serve as one more confirmation of the correctness of our 
approach. 

APPENDIX B. Power law asymptotics of equilibrium PDF 

In this Appendix we evaluate the main value of the integral 

00 

/= / ^l\k\ a+2 exp(-tkx), (B.l) 
J 2ir 

—00 

which gives the asymptotics of the equilibrium PDF, see Eq. (4.8). We present (Bl) 

as 

J = 2Re/ 1 , (B.2) 

where 

00 

/Hk 
^k a+2 exp{-tkx). (B.3) 
2ir 



In order to evaluate Ji, we pass to the complex plane. Since the integral over the 
closed contour shown in Fig. 8 is equal zero, we get 

oo 

/dk f dk 

— k a+2 exp(—ikx)=iexp(—ia7v/2) / — k a+2 exp(— kx) = 

— ico 

i exp(—iaii/2) . 



27TX a + 3 

and 



T(a + 3) (B.4) 



I= a±m%±*>. (B ,) 



7TX a+3 



With the use of Eqs. (4.8) and (B.5) we get power law asymptotics of the equilibrium 
PDF, see Eq. (4.9). 
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FIGURE CAPTIONS 



Fig.l. Equilibrium PDFs for non - linear Brownian and Cauchy oscillators in the 
linear (above) and semi - logarithmic (below) scales. 

Fig. 2. Equilibrium characteristic function for non - linear Levy oscillator in the lin- 
ear (above) and semi - logarithmic (below) scales, a = 1.5. Characteristic function 
(4.1) is evaluated by the use of Eqs. (4.1) - (4.7). 

Fig. 3. Equilibrium PDF and its asymptotics in the linear (above) and semi - log- 
arithmic (below) scales, a — 1.5. The PDF obtained by inverse Fourier transform 
of a power series (4.1) is indicated by solid lines. Dashed lines indicate power law 
asymptotics (4.9). 

Fig. 4. The profiles of equilibrium PDFs for different values of the Levy indexes. 

Fig. 5. Comparison of analytical estimates with the results of numerical simula- 
tion based on numerical integration of the Langevin equations for the Brownian 
(above) and the Levy (below) oscillators. The Levy index is 1.1. At the left: typi- 
cal trajectories, N is the number of steps. At the right the semi - logarithmic scale 
is used. The designations are as follows. Thick solid line 1 and dotted line 2 show 
the potential well and its curvature, respectively, in conventional units. Boltzmann 
distribution is indicated by solid line 3, the PDF obtained in numerical simulations 
is shown by black points, solid line 4 shows the power law asymptotics (4.9). 

Fig. 6. Comparison of analytical estimates with the results of numerical simulation 
based on numerical integration of the Langevin equations. The black points in- 
dicate the slope 7 of the power law asymptotics of equilibrium PDFs obtained in 
numerical simulations for different Levy indexes. Dotted line depicts theoretical 
value 7 = a + 3. 

Fig. 7. Comparison of analytical estimates with the results of numerical simulation 
based on numerical integration of the Langevin equations. The Levy index is 1.2. 
Thick solid line 1 and dotted line 2 show the potential well and its curvature, 
respectively. Solid line 3 indicates PDF obtained by inverse Fourier transform 
from Eq. (4.1). Black points indicate PDF obtained in numerical simulation. 

Fig. 8. Contour of integration for evaluation I 1: see Eq. (B.3). 
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